298 


N91-23007 


PRE-TERRESTRIAL ORIGIN OF "RUST" IN THE NAKHLA METEORITE. 

Susan J. Wentworth 1 and James L. Gooding 2 , ^kheed/ESC, 2400 NASA Rd. l, Houston TX 
77058. 2 SN21/Planetary Science Branch, NASA/Johnson Space Center, Houston, TX 77058. 

INTRODUCTION. As part of our campaign to identify and rigorously test evidence for aqueous 
geochemistry in shergottite, nakhlite, and chassignite (SNC) meteorites, we previously reported evidence for 
pre-terrestrial origin of calcium carbonate associated with silicate "rust" in Nakhla [1]. Here we preset 
quantitative elemental compositions and summarize textural evidence for pre-terrestrial origin of the rust. We 
denote the material in question as "rust" because its phase composition remains unknown. The rust is probably 
the same as the iddingsite-like material reported by Bunch and Reid [2] and might be related to pockets of 

material called "brownies" by Papanastassiou and Wasserburg [3]. _ , 

SAMPLES AND METHODS. Although our previous work on Nakhla [1,4,5] has also mvolved untreated 
interior and exterior chips obtained from the British Museum (Natural History), all data reported here were 
obtained from polished thin sections prepared from Nakhla specimen, USNM-426 (Smithsonian Institution). 
Scanning electron microscopy (SEM) and energy-dispersive X-ray spectrometry (EDS) foUowed our previously 
described procedures and included mineral-standardized, quantitative EDS analyses. Relative to convention 
wavelength-dispersive microprobe analyses, our method permitted beam-spot sizes of 1-10 /im, as necessitated 

by the rust occurrences. . , , . . . . ( 

RESULTS AND DISCUSSION. MICROSTRATIGRAPHY. Compelling evidence for pre-terrestnal origin ot 

the rust is found as rust veins truncated by fusion crust (Fig. 1) and preserved as faults in sutured igneous 
crystals (Fig 2). Rust veins that approach the meteorite’s fusion crust become discontinuous and exhibit vugs 
that suggest partial decrepitation; no rust veins that penetrate fusion crust have been found. Because the rust 
probably contains volatile compounds (see below), it is reasonable to expect that heating near the ablation 
surface (formed during atmospheric entry to Earth) would encourage devolatilization of the rust. Hence, 
absence of rust veins in fusion crust and vugs in rust veins near fusion crust (but not at distance from fusion 
crust) clearly imply that the rust existed in the meteorite before atmospheric entry. At least one unmistakable 
example of a vein post-dated by later endogenous events occurs as a faulted rust vein enclosed in an olivine 
crystal (Fig 2) It is obvious that the rust vein has been broken by a fault but the two disjointed sections of the 
vein are separated by clear olivine rather than an open fracture. We infer that, after faulting, the fracture was at 
least partially healed by elevated pressure. It is difficult to conceive any reasonable natural processes that 
would break and suture an igneous mineral after the meteorite arrived on Earth. Accordingly, disturbance of 
the rust vein must have occurred before the meteorite arrived on Earth and, most likely, on the Nakhla parent 
planet COMPOSITION. Our results for the average elemental composition of the rust (Table 1) gener y 
agree with those of Bunch and Reid [2] except for Al and Na. Compared with the earher data [2] we find 
much lower Al. Because our analyses of Al in Nakhla plagioclase and pyroxene agree with the respective data 
of Bunch and Reid [2], no systematic analytical error is suspected in our Al data. Our Na results for plagioclase 
agree with those of Bunch and Reid [2] but are systematically high relative to reference samples that contain 
very minor to trace Na. We conclude that the rust actually contains < 1% Na 2 Q. Low analytical totals, both m 


Table 1. Elemental compositions (weight percent) of rust in Nakhla determined by electron probe microanalysis (all iron calculated as 
Fc 20^ EDS background problem makes our Na results systematically high; nr = not reported). __ 
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Figure 2. SEM photomicrographs (secondar^.ectron images) o. 

Frame (B) is a detailed view of the area marked by an arrow ( ) 
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A MODEL FOR CRUSTAL SUBDUCTION BY LARGE IMPACTS; R. W. Wichman and P. H. Schultz. 
Dept of Geological Sciences, Brown University, Providence, R.l. 02912. 

Introduction: In the standard model of crater excavation, ballistically ejected material 
represents only about half the volume of the transient cavity (1). The other [>alf corres. ponds to 
downwardly displaced, shock-compressed material driven beneath the transient cavity (e.g., l, d, 
3, 4, 5). In large craters, the final crater then forms by the collapse of this transient .cavity ' wrth uplift 
and inward flow of the shocked, displaced material comprising the crater loor (6 ). q Although this 
idealized model of crater formation fits most well-preserved planetary impact st fu c tu re^ it imph y 
assumes an elastic halfspace beneath the target surface which may b ®' nap ^ 
laraest earlv basin-forming impacts. For these very large impacts, the depth of the transient cavty 
may excede the lithospheric thickness and, at least on Mars, such cavities apparently can interact 
with underlying viscous mantle regimes during basin formation (7). The extreme size of _these crater 
ing events also challenges some of the assumptions concerning cavity growth and collapse extrapo- 
lated from smaller structures (8). In this abstract, we propose that viscous deformation J^eathvery 
larae impacts can allow emplacement of vertically displaced crustal material in the mantle, and 
speculTteOT^he^nfplications such -impact subduction- might have i" subsequent 

Subduction Model: Although material flow fields result in ballistic trajectories for nAh* 
surface regions of the transient cavity, a full target section is preserved under the ^nter of the 
impact. In the case of an elastic half space, compression of this section against ur^fonmd 'rocks at 
greater depth enhances lateral flow, disrupting the column and spreading it across the base o the 
transient cavity. If viscous flow occurs beneath the impact, however, downward displacement of this 
crustal section is partly accomodated by lateral flow in the mantle below. This process transfers 
impact deformation from the lithospheric section to the mantle and results 

shocked crustal material into rocks of the upper mantle or asthenosphere. Although later dynamic 
rebound might limit the depths such material could reach, rebound uplifts the region i beneath the 
crater as a whole (not the displaced elements of the transient cavity alone) and the sub-impact 

crustal section initially should stay in the mantle. amn ,„ 0 

Viscous deformation of the mantle during impact is thus a necesary condition for the emplace- 
ment of crustal sections at depth. The probability of such deformation can be evaluated # by compar- 
ing the duration of the impact event to the Maxwell time (Tm) of the mantle: Tm= t/2| x e where t is 
applied shear stress, ix is shear modulus (- 10* MPa) and e is strain rate of deformation. The Maxwell 
time is defined as the time required for viscous creep under stress to equal elastic : strain (9) _, conse 
quently, viscous behavior occurs when deformation times are greater than Tm. 
tially elastic for timescales less than Tm (9). If an impact generates shear stresses of over 10 Motors 

(10 3 MPa), strain rates in an olivine mantle range from 10-* to 10 /s for ' TSfCj 1 *-.' Sdc RpSusJ 
800-1 000°C (10). These values translate to Maxwell times on the order of 1-100 seconds. Because 

mantle flow requires that the duration of impact exceed Tm, only lar ge i. „ ,oiSice 

impactors (which have impactor penetration times of several tens of seconds (8)) are likely to induce 

such beneath an impact depends on the impact angle and the 

duration of the impact relative to Tm. If we define d as the depth of the transient cavity (roughly the 
penetration depth of the projectile into the target), for near-vertical >60 fmnlrt is^hen 
proximated by the projectile diameter (D P ) but d decreases significantly as the f 

reduced to 5° (11)! While rare, near-vertical impacts are not improbable and are the most likely to 

emplace material at depth in the mantle. For a near-vertical impact with 3 fh^mnar? cavitv ex 
therefore, we expect viscous mantle deformation comparable to the size of the impact cavity ex- 
tending to depths^ of -D P beneath the base of the transient cavity. Since the thickness of the b °wn 
drivencore isof the same scale as the mantle deformation, the crustal section remains near the base 
of the transient cavity and is likely to be embedded in the basin floor after dynamic reb ° u " b ;J{ * b ® 
duration of the impact is significantly greater than Tm, however trth ° sph en c rn at © n ai can p e ^ ® 
the mantle to depths of several D P . We propose that mantle flow will engulf this d ^placed crus a 
section outright with depths of crustal burial in the range of - 0.5D P to 2D P . This burial ot crustal 
material in the mantle is reminiscent of terrestrial plate subduction and f or a pro|ec t||e 200 km in 
diameter, such ’’subduction” could bury crust to depths of between 100 and 400 km. 

The requirement of large, low-velocity impacts for this subduction mechanism tb ® ®** e " l f 
of this process in planetary history. First, impacts of sufficient size are restricted to the penod of 
basin-forming impacts before -3.7 Ga. The low impact velocity, however, is a stronger constraint on 
the occurrence of impact subduction, since only a few planetary impactor populations P ®J' I T'' 1 t 

at 5-6 km/s The planet most likely to have experienced such collisions is Mars, where impact 
%£££ rang? down t? -5 km/s for co-orbiting, heliocentric objects (12) Tbe J " e " c n e 
of these basins with a time of predicted high mantle temperatures ( 13 > a ' s ° ^° r , S c^ 
mechanism. On the Earth and Venus, impact velocities range from 15 7 40 km/s and 17-44 km/s 
respectively (8,14) and impact subduction is much less likely. Although viscous deformation beneath 
the impact is still possible at velocities of 15 km/s (for mantle temperatures on the order of 
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1300-1400°) the increased extent of of vaporization and melting beneath the impact may preclude 
preservation of a lithospheric section in the mantle during cavity collapse. Some contamination of the 
mantle by projectile or crustal components might still be expected, however. On Mercury, the 
predicted impact velocities for comets and earth-crossing asteroids (-34-44 km/s (8)) probably 
preclude impact subduction. Nevertheless, impactors from the postulated population of Vu can as- 
teroids (15) should have had much lower velocities, so impact subduction is also possible for Mer- 
cury The probability of subduction on the Moon (impact velocities down to -6 km/s (12)), depends 
more strongly on the state of the lunar interior. Thinner lithospheres associated with a magma ocean 
might permit early subduction events, but the growth of the lithosphere over time would inhibit mantle 

flow and prevent later subduction. » . . . 

Implications: Unlike the more continuous subduction of oceanic plates observed on tartri, 
impact subduction is randomly located and episodic on a global scale. Hence only random and iso- 
lated regions of a planetary mantle can be modified by this process. Long-term effects of oceanic 
subduction such as repeated passage of melts into island arcs or convection beneath back arc basins 
are thus unlikely to occur for impact subduction. In further contrast, mantle cooling associated with 
subduction of cold lithospheric plates should not occur with impact subduction due to both shock 
heating and shear deformation in the subducted section. To first order, this section is presumed to be 
in thermal equilibrium with normal mantle temperatures at the time of emplacement. Although such a 
view is over simplistic, the mantle evolution sequence presented below can be regarded as indicative 
of the relative time scales required to achieve various mantle states. .... , . , 

Impact subduction can potentially influence mantle evolution in two ways. Injection of crustal 
radiogenic elements into the mantle could affect the long-term thermal history, whereas the introduc- 
tion of crustal volatiles could affect the melting sequence. In the first case, a crustal block would 
begin to melt in -4-6 Ma for subduction to 10 kb pressure with an initial post-impact temperature of 
1000°C. Total equilibrium melting then would occur in 9-11 Ma. Significant mantle melting is unlikely 
to result in this time, since -16 Ma are needed to achieve Iherzolite melting temperatures outside a 
subducted crustal block. The melting of embedded crustal material, however, will influence the long 
term thermal evolution of the surrounding mantle. Since a fraction of melt is trapped along grain 
boundaries during porous flow, we can approximate the mantle composition after crustal melting by 
mixing a disseminated crustal component into the mantle. For mantle-crust ratios of 100:1, such 
mixing can double the abundance of heat producing elements in the mantle and these added heat 
sources eventually can induce mantle partial melts some 100-500 Ma after the subduction event. 

The subduction of volatile concentrations could produce mantle melts on much shorter time 
scales. Although water is unlikely to be a major constituent of the crustal section as a whole, water or 
ice may be concentrated in near surface regions. For a volatile-rich regolith 500 m thick with 25 /o 
porosity, subduction under a 100-km radius projectile can subduct over 1000 km 3 of water. If the 
projectile caps the subducted section and drives it into the mantle, this volatile phase ma y es " 
cape into the transient cavity and will be trapped instead near the top of the crustal section. Addition 
of such a vapor phase to surrounding mantle compositions significantly reduces the solidus tempera- 
tures and at 1000°C and 10 kb, can initiate immediate mantle melting. If carbonates are present in 
the martian regolith, the associated fluid-rich and volatile phases could possibly acheive a kimberlitic 
character 

Conclusions: Large, low velocity impacts may inject significant crustal sections into a planetary 
mantle but this process will be most efficient if the mantle yields viscously around impact-driven 
subsidence. Such behavior is most likely before 3.7 Ga on Mars, but also may have occurred on 
Mercury or the early Moon. The depth of subduction is dependent on the relative scale of impact and 
mantle flow regimes, but can achieve depths of over 100—200 km for projectiles over 100 km in 
radius. The effects of such subduction on mantle evolution are unlike those observed in terrestrial 
subduction zones and primarily reflect the effects of subducted volatile and radiogenic isotope con- 
centrations Escape of vapor into the mantle should produce kimberlite-like mantle melts soon after 
impact Crustal melts develop some 5-10 Ma after impact and enrich higher mantle regions in 
radioqenic isotopes. Finally, isotopic heating of this enriched mantle may lead to renewed mantle 
melting several hundred million years after the original impact event. Such a mantle melt sequence 
may fit the general sequence of highland volcanism observed on Mars where explosive patera vol- 
canism evidently preceded formation of most of the basaltic shields and ridged plains (16). 
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LARGE SCALE COMPRESSION STRUCTURES IN THE ERIDANIA-PHAETHONTIS REGION: MORE 

EVIDENCE FOR POLAR WANDERING. R. W. Wichman and P. H. Schultz, Dept, of Geological Sciences. 
Brown University, Providence R.l. 


INTRODUCTION: Although most tectonic structures on Mars can be associated with 
to either Tharsis formation (1,2, 3, 4) or basin modification (4,5,6), some structures in the martian 
highlands appear to be independent of both stress sources. Two ma|or systems of ^ 

rival the lobate scarps of Mercury in size and indicate that a ,' arg ® : j^'ona -scale compression^ 
event occurred in the Terra Cimmeria/Terra Sirenum region (7). These structures c annot I M ex 
Dlained bv either Tharsis-centered or basin-centered regional stress models, but may reflect defor 
mation due to polar wandering as predicted by Melosh (8). If the Tharsis Province was not so interest- 
ing tectonically, the size and isolation of these features would have ! attracted I immediate att e ntion and 
possibly stimulated discussions comparable to those of the mid-1970 s for Mercury. In abstract, 
we use V crater deformation to estimate the amount of crustal shortening across these scarp/ridge 
systems and reconsider the possible sources the causative compressional stress field. 

FEATURE DESCRIPTION: The ridge systems described here are Eridania Scopulus and a system 
of scarps and ridges located near the Copernicus impact basin The systems are of s '™^ ar ag ® 
with crater dates (derived from linear crater counts (6,9)) of -381 ± 270 and -341 ± ^240 _N(>5)/1 10 
km 2 , respectively. The Copernicus ridge system is at least 840 km in length, f 

Scopulusextends over 1000 km. Shadow measurements indicate typical scarp and ridge n h ® ^ , nH 
500-700 m with a maximum relief exceeding 1 km. The systems are thus C0 ^P a ^ de h in f' d f 
scale of deformation to the mercurian lobate scarps, which are typically several ^ u ^ gr ® d ^ 
kilometers in length and range in relief from a few hundreds of meters to one or two kilometers 
(10,11). The martian systems differ, however, in detailed morphology from ’° b f l 

scarps The mercurian scarps are typically more rounded and flatter 'in i proWe thanthe systems are 
tures; and while the mercurian scarps occur as long, individual features the martian systems 
continuous bands of deformation 100 to 200 km wide containing between 2 and 5 distinct, parallel to 
suboarallel structures. In general, these structures are scarp-like with a backslope away from the 
scaro on the uplifted side. Where they cross the floors of large craters (>20 km diameter) or inter- 
crater plains units, scarp expression is subdued and more rounded. In some plains regions the scarps 
then merqe into wrinkle ridge systems of the same orientation. Finally, the martian scarp systems 
colEckte wfth or form apparent linear topographic highs (figure 1). Such behavior is not characteristic 
of The %rcuriw tobate sois and maj reflect the combined effects of uplift across individual struc- 

tUre CR n u^AlT S^O^TEN^T jgT The nature and degree of deformation in the martian ridge/scarps are 
best revealed by the modification of craters. Deformation most commonly joins SC ^V S ™9>}er 
relief on each side of the crater with a linear, asymmetric ridge , cut ,ng ' ^ J°/J£ s ujf t differs 

floor. This pattern resembles crater deformation along the mercurian lobate scarps (1 1), but d ffers 
from the typical deformation associated with smaller wrinkle ndges, which tend to curve around^ a 
crater interior instead of cutting directly across it (4,12). A less common modification 'along I JJj® ma £ 
tian ridge/scarps preserves partial craters on the elevated side of a scarp but does P^®™® lt J® 
corresDondina down-thrown rim sections (figure 1 ) . Three craters appear foreshortened due to 
compressional deformation: one of these also shows an offset rim along the » scarp » 
mation of the Guido d' Arezzo crater by the Vostok scarp on Mercury (11 ) • For initially Cl £ c ^ ar ® rat ® f r ® 
with no extension perpendicular to compression, the difference between the long and short axes of a 
deformed crater provides a measure of the regional shortening across the structure. Such rneasure 
ments for these three craters indicate on the order of 1 .5-3 km of crustal shortening across ^indi 
scarp features. This is comparable to the values derived by Strom et al (10) for Mercury, where 
approximately 1-2 km of shortening is estimated for individual mercunan scams. 

The linear nature of deformation in crater interiors and the preservation of uplifted partial crat 
are consistent with compression by deep-seated thrust or reverse ,; au >‘'" g ^ s h or ten- 

served scarp relief helps to constrain fault models for comparison with the observed cra * er 
ina For a reverse fault with a dip of 60° and a scarp 500 m in elevation, we would expect -290 m of 
shortening, a value some 5 to 10 times smaller than that indicated by the deformed craters. AUerna- 
tivelv a thrust fault with a dip of 25° and scarp relief of 500 m would indicate -1000 m of shortening 
consistent with the estimated crater deformation. Therefore, a thrust i n Th° e d i r id 0 ania and Cooer 
appears consistent with the observed deformation and the parallel s ® ar Pf J" ' th ®. 1 ^ ' d ?hnrtTnina over 
nicus-Newton svstems might represent systems of associated thrust sheets. If the shortening 
each thrust is on the order of 1-2 km, in keeping with the crater deformation, such systems could 
accomodate up to 5 to 10 km of regional shortening perpendicular to the scarp trends 

DISCUSSION: The timing of this compressional event can be denved from the c !?j® r ag ®^ 
two systems and is broadly correlated to the time of Tharsis foimatioMfigure ^ 
error bars, however, these age determinations encompass both the time of an ? l ®5J. t J ha ^ sl ® 
and the later main phase of Tharsis radial fracturing. Although coincident in time with Tharsis centered 
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reformation the scaro/ridae systems occur over 70° from the center of Tharsis. Of the Tharsis load 
models developed by Banerdt et al (3), only the isostatic load model could provide compressive 

stresses at this P distance from Tharsis. Based on the observed def °™ at '° n 0 1Sm°th?^£ls Mrter 
Banerdt et al argue that this model only applies to regions i less rthan 40 from the Tharsis cert 
Furthermore the two scarp/ridge systems occur at a distinct angle to each other, nenner ^ iren 
consistent with the Tharsis stress orientations. Because the systems are well re ™° v ® d Jl°™ 
nearest recoanized large impact basins (Hellas, Isidis), basin-centered stresses also se ® , T ,s 
gible Even if a large buried basin did exist in the region, current models of basin-centered deform 
Son (6) could not account for the observed pattern of S "^ na , streS s field inde- 
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Figure 1. Parts of the Copernicus system (a) and Eri- 
dania Scopulus (b) ‘ , . 

scarp formation. 

• V. M 


Arrows indicate craters modified by 



pendent of both Tharsis and basin-controlled 
deformation thus appears responsible for the 
formation of these ridge/scarp systems. The 
global contraction mechanism proposed for 
Mercury (10,11) seems unlikely on Mars given 
the early timing and regional scale of the 
ridge/scarps as well as the widespread contem- 
porary extension associated with Tharsis. Polar 
wandering, however, can simultaneously 
produce regions of compressive and tensile 
stress at the several kilobar level (8), thereby 
producing the scale of deformation indicated by 
the crater shortening data. The location and 
orientation of the scarp/ridge systems are con- 
sistent with polar wandering in response to Thar- 
sis development (7), as is the correlation in age 
of ridge formation with initial Tharsis activity. 

CONCLUSIONS: Large ridge/scarp systems 
in the Terra Cimeria/Terra Sirenum region ap- 
pear to reflect a major regional event of com- 
pressional deformation with regional crustal 
shortening of up to 5-10 km. Although contem- 
porary with Tharsis fracturing, ridge system for- 
mation appears to be independent of Tharsis 
stress fields and is more consistent with defor- 
mation expected in the polar wander scenario of 
Schultz and Lutz (7). 
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